







02 449200 





NBS SPECIAL PUBLICATION /13 


U.S. DEPARTMENT OF COMMERCE /National Bureau of Standards 


Evaluating Thermal Fire 
Detection Systems 
(SI Units) 











he National Bureau of Standards’ was established by an act of Congress on March 3, 1901. The 

Bureau’s overall goal is to strengthen and advance the nation’s science and technology and facilitate 
their effective application for public benefit. To this end, the Bureau conducts research and provides: (1) a 
basis for the nation’s physical measurement system, (2) scientific and technological services for industry and 
government, (3) a technical basis for equity in trade, and (4) technical services to promote public safety. 
The Bureau’s technical work is performed by the National Measurement Laboratory, the National 
Engineering Laboratory, the Institute for Computer Sciences and Technology, and the Institute for Materials 
Science and Engineering . 


The National Measurement Laboratory 


Provides the national system of physical and chemical measurement; e Basic Standards” 
coordinates the system with measurement systems of other nations and e Radiation Research 
furnishes essential services leading to accurate and uniform physical and © Chemical Physics 
chemical measurement throughout the Nation’s scientific community, in- e Analytical Chemistry 
dustry, and commerce; provides advisory and research services to other 

Government agencies; conducts physical and chemical research; develops, 

produces, and distributes Standard Reference Materials; and provides 

calibration services. The Laboratory consists of the following centers: 


The National Engineering Laboratory 


Provides technology and technical services to the public and private sectors to © Applied Mathematics 
address national needs and to solve national problems; conducts research in e Electronics and Electrical 
engineering and applied science in support of these efforts; builds and main- Engineering” 

tains competence in the necessary disciplines required to carry out this Manufacturing Engineering 
research and technical service; develops engineering data and measurement Building Technology 
capabilities; provides engineering measurement traceability services; develops Fire Research 

test methods and proposes engineering standards and code changes; develops Chemical Engineering” 

and proposes new engineering practices; and develops and improves 

mechanisms to transfer results of its research to the ultimate user. The 

Laboratory consists of the following centers: 


The Institute for Computer Sciences and Technology 


Conducts research and provides scientific and technical services to aid 


Programming Science and 


Federal agencies in the selection, acquisition, application, and use of com- Technology 
puter technology to improve effectiveness and economy in Government © Computer Systems 
Operations in accordance with Public Law 89-306 (40 U.S.C. 759), relevant Engineering 


Executive Orders, and other directives; carries out this mission by managing 
the Federal Information Processing Standards Program, developing Federal 
ADP standards guidelines, and managing Federal participation in ADP 
voluntary standardization activities; provides scientific and technological ad- 
visory services and assistance to Federal agencies; and provides the technical 
foundation for computer-related policies of the Federal Government. The In- 
stitute consists of the following centers: 


The Institute for Materials Science and Engineering 


Conducts research and provides measurements, data, standards, reference © Ceramics 

materials, quantitative understanding and other technical information funda- e Fracture and Deformation 2 
mental to the processing, structure, properties and performance of materials; e Polymers 

addresses the scientific basis for new advanced materials technologies; plans e Metallurgy 

research around cross-country scientific themes such as nondestructive e Reactor Radiation 
evaluation and phase diagram development; oversees Bureau-wide technical 

programs in nuclear reactor radiation research and nondestructive evalua- 

tion; and broadly disseminates generic technical information resulting from 

its programs. The Institute consists of the following Divisions: 





"Headquarters and Laboratories at Gaithersburg, MD, unless otherwise noted; mailing address 
Gaithersburg, MD 20899. 

*Some divisions within the center are located at Boulder, CO 80303. 

‘Located at Boulder, CO, with some elements at Gaithersburg, MD. 








| 


Evaluating Thermal Fire 
Detection Systems 
(SI Units) 





David W. Stroup 
David D. Evans 
Phyllis Martin 


National Bureau of Standards 
National Engineering Laboratory 
Center for Fire Research 
Gaithersburg, Maryland 20899 


U.S. DEPARTMENT OF COMMERCE, Malcolm Baldrige, Secretary 
NATIONAL BUREAU OF STANDARDS, Ernest Ambler, Director 


Issued April 1986 


Library of Congress Catalog Card Number: 86-600520 


National Bureau of Standards Special Publication 713 
Nati. Bur. Stand. (U.S.), Spec. Publ. 713, 557 pages (Apr. 1986) 
CODEN: XNBSAV 


U.S. GOVERNMENT PRINTING OFFICE 
WASHINGTON: 1986 


For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, DC 20402 














TABLE OF CONTENTS 


eS OMMLAB ILE Si rewolchlelerodslethshoiadolcNcleVolonelclotel ofelclctet) evellolte eRelclcdoielohotele 5O0 0000 
LTS TOK PGURESlereteneleloielolel ofekonelcletehokollslclelenslslevehchel> 560000000000 S5000005C 
NOMENGIEAT URE evetertetaralenetenalcleisl ichehoiaistelaleiclekclelciielelelcllcreliaiatelonele/elaielicliele +) 500 
/NSEIETCEYSIE 9 05000000000000000000006000560606050000000000000000000000 
Iho  TUSHEO)DUGIILO) 6 Go 66 00G0uD0000000 66000000 900000000000000000000000006 
PRODESCREPTION OF THE MODEL. 2 oleic 100s ooWltsetle at IF DR PEATE 6 658 
oil Grerserearedtoete WilssBo0900600000000000000000000000050000000000000 
2.2 Temperature and Velocity Correlations............ccccccceeees 
2.3 Detector Response Model........ccccccvcescecvees 5000000000000 
2.4 Velocity Correlation for r/H < 0.3..... pedis: slave’ Wevereveioleveie eaves 
3. IMPLEMENTATION OF THE MODEL: COMPUTER PROGRAMS.........cccccceees 
3.1 General Program Structure.......cccecccecrcccccsceccvccceces O00 
Bo Weal ABORT. 9 5000000600000000000000000000000000000000000 S000 
3.2.1 Data Input Ser eton eeeeoee?es eeeeoeeeeseeeeeeeeeeeeeeeeeeeee @ 
3.2.2 Calculation Initialization........ 5000000000000 90000000 
3.2.3 Conversion and Output Section..............0e. 500000060 
3.3 Activation Time Equation Solution: 
Ipsiscyel TEM ENSENEOEFALG 6 6 G0600000000000000000000000000000000000 
3.4 Activation Time Equation Solution: 
Rate of Temperature Rise............. 900000000000000000000 
4. IMPLEMENTATION OF THE MODEL: TABLES OF RESULTS.......... 500000000 
5. EXAMPLES @eeoeoeeevee@ @eeeeeoeeeee eee 88 8 6 @ eee ee e@eeeoeoee2e e e ee @eeeoe ee 
Soll Weta he WV WASs ooodandb 000000 UODOO0UO OO ODO UOOG og0d0cCC00N . 
5.1.1 Example 1: Comparison with NFPA 72E 
PablewReSuiteSireiatalerelelelaalcialelelelsl ele Sogo DOD OD DOO DDOOONOCO 
5.1.2 Example 2: Fixed Temperature vs. 
Rate of Rise.........eeeeee 600000 5000000000000000000 
5.2 Example 3: Using the Computer Program...........ceesoes 900000 
6. SUMMARY @eeeoeeees eee @ @eeoeoeoeveeveeveeeee88 888 8 &@ & @ e @eeeoeaeeeoeeeees eee eee ee e e 
7. REFERENCES........... pealatateteevatetatelevalatencterelel sires Sa060000N S60 00000 : 
APPENDIX A - DETACT-T2 (FORTRAN VERSION).......... SreleKeiekeleialeiels)ereKeie 


iii 


TABLE OF CONTENTS (continued) 


Page 
APPENDIX B - DETACT-T2 (PC BASIC VERSION)............0-. HelleloteHone! tent a-OS 
APPENDIX C - DETECTOR ACTIVATION TIME TABLES: 
FIXED TEMPERATURE....... ScoCDD0GCODODCOCOOOOO. DIL 
APPENDIX D - DETECTOR ACTIVATION TIME TABLES: 
RATE OF TEMPERATURE RISE...........+e0e2ee00-- JOL 
NOTE: Appendixes C and D are available in English units in another 


version of this report. 


iv 


LIST OF TABLES 


Page 
Table 1. Heat Release Rates of Warehouse Materials.................. 23 
Table 2. Maximum Heat Release Rates from 
Fire Detection Institute Analysis...........cccccccccscecee 27 
Table 3. RTI Values for any Listed Detector 
(BixedsTemperatures Type) ccccccsc ssw sede cs vcs svcesveeeesees 28 
Table 4. RTI Values for any Listed Detector 
GRatesofuRise) Type)isrs cs els cc +. s oe okeus Guheree: SR io sete Bayerore Syerenee's 29 
LIST OF FIGURES 
Page 
Figure 1. Plot of t? Fire Growth Curves............. Wevstigveuroy susie apes atots 30 
Figure 2. Plot of Dimensionless Velocity Correlation................ 31 


RTI 


AT 


NOMENCLATURE 


g/(C, T,, p,,) 

specific heat capacity of ambient air 
constant used in equations 9 and 11 
constant used in equations 10 and 11 
acceleration of gravity 

vertical distance from fuel to ceiling 
fire energy release rate 

dimensionless fire energy release rate 


radial distance from fire axis to the detector 


response time index, the product of the detector thermal time constant 
and the square root of the gas velocity used in the test to measure the 
time constant [8]. 


time 


[5 q t/S yt/5 j 


dimensionless time, t/{a 
dimensionless time for time delay for gas front travel. 
ambient temperature 

gas temperature 

temperature of detector sensing elements 


1h Ook 


(+) 


vi 


AT 


Nh 


NOMENCLATURE (continued) 


NS (T./g) a 


dimensionless temperature difference, AT/[A 

gas velocity at the detector location 

maximum gas velocity 

: 5 1/5 

dimensionless gas velocity, U/[A a H] 
P P 2 Q 0h 

proportionality constant for t™ - fire growth = Q/t 


ambient air density 


detector thermal time constant 


vii 


2/5 ,-3/5 


] 



































* 
| 
oF 
i 
ie 
{ Varn fee 
} ae 
Fe if ; 
Noe. 
ij 1 : P 
; Pl i ie 
ay yt 
ti 
Tih! 
i 
Tal 
i Pe 
i i a, 
{ i i 
i) Vt A 
Vp Da ll ten fy ; i oa Wa _ 
i / ‘ hed h \ i 

yu y y i 

Thee 

7 fat) 

+ \ : ye 
1 yu 1 





Evaluating Thermal Fire Detection Systems [S.I. Units] 


David W. Stroup 
David D. Evans 
Phyllis Martin 


ABSTRACT 


This report presents a methodology for evaluating heat detection systems 
installed in buildings. Previous work for use primarily in designing new 
thermal fire detection systems was used as a starting point. The previous work 
was enhanced and supplemented to make it more useful for evaluating existing 
systems. The resulting equations were programmed into a user interactive 
computer program. This program is available in both BASIC and FORTRAN and will 
run on mainframes as well as personal computers. In addition, a modified 
version of the FORTRAN program was used to develop an extensive set of tables 
listing detector activation times for given building geometries, detector 
characteristics, and fire growth rates. These tables are useful for quick 
evaluation of alternative heat detector installations. Finally, practical 
examples are included to illustrate the use of the tables and computer 
programs. 


Keywords: fire alarm systems; fire detection; fire detection systems; fire 
hazard assessment; fire protection; fire suppression; heat detectors; sprinkler 
systems. 


1. INTRODUCTION 


Studies of the response of heat detectors to fire driven flows under 
confined ceilings have been conducted since the early 1970's [1,2,3,4,5]!. 
Results of these largely experimental studies have been used to develop 
correlations of data that are useful under a broad range of fire conditions and 
building geometries. These correlations have been used to construct 
engineering methods to determine heat detector spacing, sprinkler response 
time, and smoke detector alarm times for industrial buildings where large 
undivided ceilings over storage and manufacturing facilities are common. The 
National Fire Protection Association (NFPA) has adopted an alternate design 
method for determining heat detector spacing which is based on some of these 


correlations. 


This alternate design method is contained in Appendix C of the NFPA 72E 
standard [6], "Guide for Automatic Fire Detector Spacing," which consists of 
tables listing maximum allowable detector spacings for given room sizes and 
fire growth rates calculated from correlations of experimental measurements. 
The tabular form of presentation is useful for specifying design requirements 
in new installations. However, the tables are inconvenient to use in 
evaluating changes to existing detector systems where the spacing is fixed. 
There appears to be a need for a method which determines detector activation 
times for specified building geometries, detector characteristics, detector 


spacings, and fire growth rates. This method would be useful in evaluating the 


1 Numbers in brackets refer to literature references at the end of this report. 


Z 


expected fire detection time for existing and proposed systems for use in 
determining available safe egress time, and the relative effectiveness of 


different detectors used at the same spacing. 


This paper presents the results of work to develop an alternative method 
appropriate to existing installations. As part of this study, the basis for 
the calculation method published in Appendix C of the NFPA 72E standard was 
determined. Alternate correlations of the experimental data, used to develop 
the results presented in Appendix C of NFPA 72E, have been used to construct a 
FORTRAN program (DETACT-T2). This computer program evaluates the response time 
of heat detector systems [7]. In this work, a modified version of the 
DETACT-T2 program was used to generate tables of detector response times for a 
wide range of building geometries, detector characteristics, and fire growth 
rates. These calculated values for response time agree to within 5 percent of 


those published in the tables contained in Appendix C of the NFPA 72E standard. 


2. DESCRIPTION OF THE MODEL 


The problem of heat detector response is solved by predicting the time 
dependent temperature of the detector sensing element up to the point when it 
is heated to the specified alarm conditions. In order to predict the 
temperature of a heat detector element, it is necessary to describe the fire 
and the fire generated environment to which the detector is exposed. The 
primary fire related parameters which influence the response of heat detectors 


are the gas temperature and velocity. Once the fire-generated environmental 


characteristics are known, a model of the response of the thermal sensing 
element in the detector is employed to calculate the detector element 


temperature as a function of time. 


2.1 Characteristic Fire 


Extensive full-scale tests have been conducted to determine the response 
of heat detectors [1-4]. These tests have studied two classes of fires. The 
first of these are steady fires. These fires have energy release rate 
histories which do not change with time. Results from work using this type of 
fire may be used to calculate heat detector response if the quasi-steady state 
assumption is invoked. In this case, a varying energy release rate history is 
approximated as a series of steady fires. In this type of analysis, the 
transport time of the gases to the detector is ignored. This usually results 
in a predicted actuation time which is earlier than that encountered in actual 


practice. This may or may not be acceptable depending on the case. 


Another class of fires for which research results are available are called 
"t-squared" (t?) fires. These fires have energy release rate histories which 
increase proportionally with the second power of time from ignition. This 
class of fires was used to develop the tables in the NFPA 72E document and will 


be used in this paper. 


Appendix C of the NFPA 72E standard lists three fires which have t? type 


fire growth behavior. A slow developing fire is defined as one which would 


take 600 seconds (10 minutes) from the time of flaming ignition until it 
reaches a heat release rate of 1055 kW (1000 BTU/s). A medium growth rate fire 
reaches a heat release rate of 1055 kW (1000 BIU/s) approximately 300 seconds 
(5 minutes) after flaming ignition occurs. The fast developing fire requires 
150 seconds (2.5 minutes) to reach a heat release rate of 1055 kW (1000 

BTU/s). In addition, a fourth t? fire has been developed in this report. This 
fire has been defined as one which takes 75 seconds (1.25 minutes) to reach a 
heat release rate of 1055 kW (1000 BTU/s). These four fire growth curves are 


plotted in Figure 1. 


Table 1 is a reproduction of a table in NFPA 72E. It lists various 
warehouse commodities and their maximum fire heat release rates per unit floor 
area covered by the commodity. The time required for a fire, occurring in a 
particular commodity, to reach a heat release rate of 1055 kW is also listed. 
If possible, each commodity is identified with a specific fire growth behavior 
- slow, medium, fast, or ultrafast fire. Table 2 is a reproduction of another 
table found in Appendix C of NFPA 72E. It lists the maximum heat release rates 
for various furnishing type fuel items. Both of these tables are useful for 


estimating fire growth rates that may be expected from various fuel packages. 


2.2 Temperature and Velocity Correlations 


Values of the time dependent gas velocity and temperature for t* fires may 


be obtained from work by Heskestad and Delichatsios [4]. From experiments and 


use of dimensional arguments, they developed correlations of maximum 


temperature rise and velocity for unconfined ceiling layer flows as a function 
of radial distance from the axis of fire plume impingement. These correlations 


of dimensionless maximum temperature rise and dimensionless maximum velocity 


ares: 
0, 35 = (). 
aT, = (1) 
3S - 0.954 (1 + r/H) We 
0.188 + 0.313 r/H P oa > GD. 
v,/ aT, = 0459 (c/Hio (2) 
where 
oa Ss  // falls q ifs wo! 5; 
us oni tal/s g/5 t/5; 
aT, er taz/5 (1_/g) a2!5 43/5) 
A = g/ le Ty, Poo] 
a = Q / if 
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(to), = 0.954 (1 + r/H) 


These equations may be used to predict the fire-generated environment (gas 


temperature and velocity) near a detector. 


2.3 Detector Response Model 


Characterization of the thermal response of heat detector and sprinkler 
thermal sensing elements is discussed by Heskestad and Smith [8], and Evans 
[9]. The model for the detector sensing element temperature is based on a 
convective heat transfer process. The first order differential equation that 


describes the rate of temperature increase of the sensing element is: 


= —— ([(T-T]) (3) 


This equation predicts the temperature change of a detector with known thermal 
characteristics to a given time dependent gas flow temperature and velocity. 
In particular, the rate of change of the detector element temperature at any 
time is equal to the gas velocity to the one-half power multiplied by the 
difference between the temperature of the hot fire generated gas flow near the 
detector and the temperature of the detector element. All divided by the 


detector RII or response time index. 


The RTI is a measure of the detector's thermal lag and can be determined 
for sprinkler or heat detector sensing elements using the plunge test [8]. 
NFPA 72E reports detector sensitivity in terms of time constants, t. At any 
given time, the time constant for a given type of detector may be calculated 


from the RTI using: 


c= Rt / vil? (4) 


Tables 3 and 4 relate heat detector listed spacings to RTI values for fixed 
temperature detectors and rate of rise detectors respectively. These Tables 


were derived from results published in references 3 and 5. 


Beyler [10] used the correlations of Heskestad and Delichatsios [4] to 
derive a new form of equation (3). He was able to solve this new form of the 
heat detector response equation analytically. Using the condition that the 
initial detector element temperature is equal to ambient temperature (T, = T), 
the equations for calculating the response of fixed temperature and rate of 


temperature rise detectors are, from Beyler: 





* *1/4 
dT, bie (4/3) (aT/aT,) (aT,) -y 
= SS eee OL =e: (5) 
dt (t/t,) (0.188 + 0.313 r/H) 
T - (0) = (AT/ AT.) Atom (ote @sa7) (6) 
s s = 2. 2 Y 
where 
* \ 1/2 * 
1/2 U AT 
Su 2 2 t 
Y=5(-> == sas — | (0.188 + 0.313r/H) 
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Using equations (1) and (2), equations (5) and (6) may be solved to predict the 
response of sprinkler and heat detector sensing elements to fire-driven ceiling 


jet flows. 


2.4 Velocity Correlation for r/H < 0.3 


While this exact mei ientte equation (3) is applicable to any r/H value, 
the dimensionless velocity correlation, previously discussed, is only valid for 
r/H values greater than 0.3. In an effort to extend the usefulness of the 
detector response calculations, an analysis was made of the available data for 
ceiling jets at r/H values less than 0.3. Results of this analysis indicated a 
lack of much quantitative information for this region. Therefore, an analysis 
to develop a rational estimation of velocities in the ceiling jet close to the 
plume axis was undertaken by studying a simple extension of the velocities in 


the plume. 


The dimensionless velocity and temperature change relations proposed by 
Heskestad and Delichatsios [4] together with equations proposed by Zukoski 
et. al. [11] to describe fire plume maximum temperature and velocity 


distributions were taken as a starting point. 


U bel ay cal/5 at! yi! 5; (7) 
AT; = AT / raz/5 (T/g) a2! H 3/5) (8) 


o -% 2/3 
Aw fei chan") (9) 


gil? y/2 (Q*)1/3 (10) 


U = 
m Oy 
These equations were used to obtain an expression of dimensionless maximum gas 


velocity in the same form as equation (2). 


vy fit, = GG au) 
Using values suggested by Zukoski et. al. [11], the dimensionless gas velocity 
is found to be a constant at r/H values of less than 0.3. The complete 
dimensionless gas velocity equation for all values of r/H is: 


: / JaT 
Ui Us) 


" 3.87 / 9.115, r/H $ 0.3 


-0.63 ee 


0.59 (r/H) r/H > 0.3 


This equation is plotted in Figure 2. The two curves merge nicely. This 
equation together with the previously discussed correlations and analytic 
solution may be used to calculate any aspect of the response of thermal 


detectors to t? fires. 


3. IMPLEMENTATION OF THE MODEL: COMPUTER PROGRAMS 
3.1 General Program Structure 


The equations discussed in Section 2 were programmed into a user 
interactive computer code called DETACT-T2. To evaluate the response of a 
detector, the user enters values for ambient air temperature, detector 
activation temperature and rate of temperature rise, detector RTI, fuel to 
ceiling distance, detector spacing, and a fire growth rate constant, a (for t? 
fires). The program solves the appropriate equation to obtain the response 
time for a fixed temperature detector and a rate of rise detector with the 
given characteristics. The fire energy release rate at each time step is also 
calculated. | 
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The computer program is available in two versions. One version, contained 
in Appendix A, is written in FORTRAN 77. The program contains approximately 
400 lines of source code and requires about 64k of computer memory in compiled 
form for execution. The source code may be compiled on any computer supporting 
an ANSI standard FORTRAN 77 compiler. The program has been compiled and run 
successfully on IBM2 and compatible personal computers operating under 
MS-pos3. A specially modified version of this program was used to generate the 


tables discussed in section 4. 


The other version of the program, listed in Appendix B, is written in 
BASIC. It consists of approximately 350 lines of source code. This program 
has also been successfully tested on IBM and compatible micro-computers. Both 
programs have essentially the same structure and use the same variable names. 


To insure convergence, all variables are in double precision. 


The program consists of three major sections. The main program reads and 
lists the required input data, performs the calculations necessary to set up 
the equations, calls the solution subroutines, and writes out the results. The 
second section is a subroutine which uses a Newton-Raphson technique to solve 
the equation describing the activation time of a fixed temperature detector. 
The last section is a subroutine which uses a bisection method to solve the 

2 IBM is a trademark of International Business Machines Corporation and 
does not represent an endorsement by the National Bureau of Standards. 


3 MS and Microsoft are trademarks of Microsoft Corporation and do not 
represent an endorsement by the National Bureau of Standards. 
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rate of temperature rise detector activation time equation. Each of these 


sections is described in more detail below. 


3.2 Main Program 


3.2.1 Data Input Section 


This section of the main program requests information from the user 
concerning the building geometry, the detector characteristics, and the fire. 
growth rate. The user may also elect to enter data in either English or 


S.I. units. 


3.2.2 Calculation Initialization 


Several preliminary calculations must be made before the activation time 
equations may be solved. These calculations are done in the second section of 
the main program. In this section, the dimensionless time, velocity, and 
temperature are calculated. Based on the user-specified detector spacing, the 
maximum radial distance of the detector from the fire is calculated. After 
these calculations have been completed, the two activation time calculation 
subroutines are called. The results from these subroutines are converted into 


activation times, and the heat release rates are calculated. 


12 


3.2.3 Conversion and Output Section 


The last section of the main program writes out the solutions to the 


equations. The results are printed in both English and S.I. units. 


3.3 Activation Time Equation Solution: Fixed Temperature 


The equation describing the activation time of a fixed temperature 
detector cannot be solved explicitly for activation time. Therefore, a 
numerical technique must be used. A Newton-Raphson technique was selected for 
use in solving the fixed temperature detector equation. This technique is 


discussed in most numerical methods textbooks [12]. 


The Newton-Raphson (or simply Newton's) method is one of the most powerful 
and well-known numerical methods for finding a root (solution) of f(x) = 0. To 
use this routine, the activation time equation is first adjusted to fit the 
appropriate form. This method converges very rapidly for the cases considered 


here. 


3.4 Activation Time Equation Solution: Rate of Temperature Rise 


A bisection technique was selected for use in solving the rate of rise 


activation time equation. This technique is also discussed in most numerical 


methods textbooks [12]. It was chosen for use in this instance because the 
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Newton method would have introduced stability problems into the calculation. 


This method also converges very rapidly for the cases considered here. 


4. IMPLEMENTATION OF THE MODEL: TABLES OF RESULTS 


The computer programs calculate detector activation times one case at a 
time. In many instances, it is necessary to evaluate a number of 
alternatives. It would be tedious and inconvenient to run a computer program 
to examine each case of interest. In light of this fact, the FORTRAN version 
of the computer program was used to generate an extensive set of tables. These 
tables list detector activation time as a function of building geometry, 


detector characteristics, and fire growth rate constants. 


Each entry in the tables was checked by first calculating the response 
time of the given detector numerically. Then the calculated time is 
substituted into the exact form of the equation, and it is solved for the 
temperature of the detector sensing element at that time. In each case, the 
calculated temperature agreed with the originally assumed one to within an 


accuracy of 0.005 percent. 


The tables cover a wide range of parameter variations. Tables have been 
developed for each of the four fire growth rates that were described earlier. 
Tables are available in both English and S.I. units. The S. I. units tables 
are contained in Appendixes C and D of this report. The English units tables 


are available in the English version [13]. 


14 


For use with S.I. units, tables cover detector spacings of 0 to 15 meters 
and ceiling heights of 1 to 24 meters. RTI values range from 0 to 400 (m-s)?/? 
for fixed temperature detectors and from 50 to 1000 (m-s)?”? for rate of rise 
detectors. Temperature differences for fixed temperature detectors range from 
5 to 80 degrees C, and available rate of temperature rise values are 8, 10, and 


14 degrees C/minute. 


The tables of detector activation times are organized according to: type 
of detector (fixed temperature or rate of temperature rise), detector spacing, 
and type of fire (ultrafast, fast, medium, or slow). Appendix C contains the 
tables listing activation times for fixed temperature heat detectors, and 
Appendix D contains the tables listing activation times for rate of temperature 
rise detectors. In each appendix, the tables are arranged in order of 
increasing detector spacing (i.e. 0 m, .5 m, 1 m, etc.). For each spacing, the 
tables are arranged according to fire type. Ultrafast is first, followed by 


fast, medium, and slow. 


The asterisks next to some of the table entries indicate that the detector 
location places it within the r/H < 0.3 region. This region is where the 
theoretical dimensionless velocity correlation, discussed in section 2.4, is 
applicable. These values are flagged because the calculation of gas velocity 
in this region is based primarily on a theoretical analysis and not on a 
combination of theory and experiment. In addition, this region is not 


considered in the NFPA 72E Appendix C. 
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Finally, just as in Appendix C of NFPA 72E, for the tables in this report, 
the fire plume is assumed to impinge on the ceiling at the center of the square 
array of detectors spaced as indicated. So for detectors spaced at 5 m (16.4 
ft), the distance of each detector from the gas impingement point would be 3.53 
m (5 42/2). Each table indicates this distance in parenthesis after the table 


spacing value. 


5. EXAMPLES 
5.1 Using the Tables 
5.1.1 Example 1: Comparison with NFPA 72E Table Results 


The first example will show that these new tables agree quite well with 
the results from the tables in Appendix C of NFPA 72E even though each set of 
tables were generated using different correlations of experimental data. 
Assume a system of detectors is installed using a 3.5 meter spacing in a room 
with a 4 meter ceiling height. A medium fire growth rate is expected. The 
detectors have an RTI rating of 150 (m-s)!/%? which corresponds to a 121.5 
second time constant. The ambient temperature is 21 degrees C and the rated 
activation temperature of the detector is 81 degrees C. This gives a 
temperature difference of 60 degrees. 

AT =T 


= dt 81 - 21 = 60 degrees C 


amb * 
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Looking in Appendix C, a detection system installed using a 3.5 meter 
spacing would detect a medium growth rate fire 5.52 minutes after ignition. If 
we next look for similar conditions in the 72E tables [NFPA, p. 72E-36], we 
find that detectors with a time constant of 100 s installed using a 3.5 m (11.5 
ft) spacing in a room with a 3.66 m (12 ft) ceiling height will activate when a 
medium growth rate fire reaches a heat release rate of 1055 kW. A medium 
growth rate fire will reach this heat release rate 300 seconds or five minutes 


after ignition. 


The detector listed in the NFPA table has a time constant of 100 s, and 
its activation temperature difference is about 55.5 degrees. These are 
slightly less than the 121.5 s and 60 degrees C, respectively, for the detector 
listed in the NBS table. Therefore, it is expected that the detector listed in 
the NBS table would have an activation time slightly higher than the one in the 
NFPA table. After extensive cross checks between the new NBS tables and the 


72E tables, the results were found to agree to within 5 percent. 


5.1.2 Example 2: Fixed Temperature vs. Rate of Rise 


This example will show how alternative heat detector installations may be 
compared using the tabular results. In particular, the differences between 
fixed temperature and rate of rise detectors will be compared. Suppose a 
detection system is installed in a room with a ceiling height of 12 meters and 
detection of a medium fire is desired. The detectors are spaced 7 meters 


apart. Currently, fixed temperature detectors with an RTI of 150 (m-s)}4? and 


7 


an activation temperature rating of 81 degrees C are installed in the room. 
Ambient temperature is 21 degrees C. (Temperature differential of 60 

degrees.) Will the use of rate of rise detectors with the same RTI as the 
fixed temperature detectors and an activation rate of rise rating of 8 degrees 
C/min yield a lesser time to activate a detector in response to a medium growth 


rate fire? 


From Appendix C, the activation time of the fixed temperature detectors, 
installed using a 7.0 meter spacing and exposed to a medium growth rate fire, 
is 13.64 minutes. The time required to activate the rate of rise detectors is 
obtained from Appendix D. The activation time for these detectors is 15.82 
minutes. From Figure 1, the heat release rate at activation of the rate of 
rise detectors would be 10,561 kW compared to 7851 kW for the fixed temperature 


detectors. 


It would take about three minutes longer to activate the rate of rise 
detectors, and the fire heat release rate would have increased by about 2710 
kW. This is interesting in light of the generally held belief that rate of 
rise detectors activate before fixed temperature detectors. This is not true 
when the detectors are located in areas with high ceilings and the fire does 
not grow rapidly. This example illustrates how quick and sometimes surprising 
comparisons may be made using the tabular form of the detector activation time 


results. 
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5.2 Example 3: Using the Computer Program 


The last example will be used to show how the computer programs work. 

Both the BASIC version and the FORTRAN version ask similar questions and 
operate in a similar manner so this discussion is applicable to both programs. 
The problem will be to calculate the activation time for fixed temperature and 
rate of rise detectors installed, using a 3 meter spacing, in an area with a 
ceiling height of 4 meters. The detectors have an RTI of 370.3 (m-s)!4%?. The 
detector activation temperature is 54.4 degrees C, and the activation rate of 
rise is 8.33 degrees C/min. Ambient temperature is 21 degrees C. (User input 
is preceded by a > and, where possible, typed in lower case.) 
ENTER 1 FOR ENGLISH UNIT INPUT 

2 FOR METRIC UNIT INPUT 
>2 


ENTER THE AMBIENT TEMPERATURE IN DEGREES C. 
>21 


ENTER THE DETECTOR RESPONSE TIME INDEX (RTI) IN (M-SEC)**1/2. 
>370.3 


ENTER THE DETECTOR ACTIVATION TEMPERATURE IN DEGREES C. 
>54.4 


ENTER THE DETECTOR RATE OF RISE IN DEGREES C/MIN. 
>8.33 


ENTER THE CEILING HEIGHT IN METERS. 
>4 


ENTER THE DETECTOR SPACING IN METERS. 
>3 


ENTER: S FOR SLOW FIRE GROWTH RATE, 
M FOR MEDIUM FIRE GROWTH RATE, 
F FOR FAST FIRE GROWTH RATE, 
U FOR ULTRAFAST FIRE GROWTH RATE, OR 
O FOR OTHER. 
>m 
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RESULTS: 


CEILING HEIGHT 


= 4.00 METERS ( 13.12 FEET) 
DETECTOR SPACING = 3 


.00 METERS ( 9.84 FEET) 
DETECTOR RTI = 370.3 (M-SEC)**1/2 ( 670.7 (FT-SEC)**1/2) 
FIRE GROWTH RATE CONSTANT = .1172E+002 JOULES/(SEC**3) 
(.1111E-001 BTU/SEC**3) 
FOR TEMPERATURE ACTUATED DETECTOR: 
ACTIVATION TEMPERATURE = 54.4 DEGREES C (129.9 DEGREES F) 
TIME OF ACTIVATION = 5.09 MINUTES 


HEAT RELEASE RATE = .1094E+004 KILOJOULES/SEC 
(.1037E+004 BTU/SEC) 


FOR RATE OF RISE ACTUATED DETECTOR: 


ACTIVATION RATE OF RISE = 8.33 DEGREES C/MIN 
(14.99 DEGREES F/MIN) 


TIME OF ACTIVATION = 3.15 MINUTES 


HEAT RELEASE RATE = .4199E+003 KILOJOULES/SEC 
(.3980E+003 BTU/SEC) 


The results show that the fixed temperature heat detector would activate 
approximately 5 minutes after the fire reaches a flaming state. The heat 
release rate at this time would be 1094 kW. A rate of rise detector would 
activate at about 3 minutes with a corresponding heat release rate of 419.9 
kW. (If English units had been selected, the input requests would have called 


for data in English units instead of metric units.) 
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6. SUMMARY 


This report presented a methodology for evaluating heat detection systems 
installed in existing buildings. Previous work for use primarily in designing 
new systems was enhanced and supplemented to broaden its application. The 
resulting equations were programmed into a user interactive computer program. 
This program is available in both BASIC and FORTRAN and will run on mainframes 
as well as personal computers. In addition, a modified version of the FORTRAN 
program was used to develop an extensive set of tables listing detector 
activation times for given building geometries, detector characteristics, and 
fire growth rates. These tables are useful for quick evaluation of alternative 
heat detector installations among other uses. Practical examples were utilized 


to illustrate the use of the tables and computer programs. 
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TABLE 1. WAREHOUSE MATERIALS” 


Time to 
Reach 
1055 kW 
(sec) 


. Wood pallets, stack 


1.5 ft high (6-122 
moisture) 150-310 


.- Wood pallets, stack 


5 ft high (6-122 
moisture) 90-190 


- Wood pallets, stack 


10 ft high (6-12% 
moisture) 80-110 


- Wood pallets, stack 


16 ft high (6-122 
moisture) 75-105 


. Mail bags, filled, 


stored 5 ft high 190 


- Cartons, compart- 


mented, stacked 15 
ft high 60 


. Paper, vertical rolls, 


stacked 20 ft high 15-28 


. Cotton (also PE, 


PE/Cot, Acrylic/ 

Nylon/PE), gar- 

ments in 12 ft high 

rack 20-42 


Cartons on pallets 
rack storage, 15-30 
ft high 40-280 


Maximum 
Heat Classification 
Release (s - slow) 
per Unit (m - medium) 
Floor Area (£ - fast) 
(kW/m?) (u - ultrafast) 





1249 m-f 


3748 f 


6814 f 


10221 f 


397 f 


Adapted from Table C-2-2.2.1(a) of Ref. 5 
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10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


TABLE 1. WAREHOUSE MATERIALS (continued) 


Paper products, 
densely packed in 
cartons, rack 
storage, 20 ft high 


PE letter trays, 
filled, stacked 5 ft 
high on cart 


PE trash barrels in 
cartons stacked 15 ft 
high 


FRP shower stalls in 
cartons, stacked 15 
ft high 


PE bottles packed in 
Item 6 


PE bottles in car- 
tons, stacked 15 ft 
high 


PE pallets, stack 3 ft 
high 


PE pallets, stack 6-8 
ft high 


PU mattress, single, 
horizontal 


PF insulation board, 
rigid foam, stacked 
15 ft high 


PS jars packed in 
Item 6 


Time to 
Reach 

1055 kW 
(sec) 


470 


190 


55 


85 


85 


75 


130 


30-55 


55 
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Maximum 
Heat 
Release 
per Unit 
Floor Area 


(kW/m?) 


8517 


2839 


1249 


6246 


1931 


Classification 
(s - slow) 
(m - medium) 
(£ - fast) 
(u - ultrafast) 


21. 


Dore 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


30. 


31. 


TABLE 1. WAREHOUSE MATERIALS (continued) 


PS tubs nested in 
cartons, stacked 14 
ft high 


PS toy parts in car- 
tons, stacked 15 ft 
high 


PS insulation board, 
rigid, stacked 14 ft 
high 


PVC bottles packed 
in Item 6 


PP tubs packed in 
Item 6 


PP and PE film in 
rolls, stacked 14 ft 
high 


Distilled spirits in 


barrels, stacked 20 ft 


high 

Methyl alcohol 
Gasoline 
Kerosine 


Diesel Oil 


Abbreviations: 


Time to 
Reach 

1055 kW 
(sec) 


105 


110 


10 


40 


FRP - Fire Retardant Plastic 


= Polyethylene 


= Phenolic Resins 


- Polypropylene 
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Maximum 
Heat 
Release 
per Unit 
Floor Area 
(kW/m? ) 


5110 


2044 


3293 


3407 


4429 


3975 


2271 


2044 


Classification 
(s - slow) 
(m - medium) 
(£ - fast) 

(u - ultrafast) 





Polysytrene 
Polyurethane 
Polyvinyl Chloride 
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TABLE 2. MAXIMUM HEAT RELEASE RATES FROM 


FIRE DETECTION INSTITUTE ANALYSIS 


Approximate 
Values 
kW 
1. Medium wastebasket with milk cartons 105 
2. Large barrel with milk cartons 148 
3. Upholstered chair with polyurethane foam 369 
4. Latex foam mattress (heat at room door) 1266 
5. Furnished living room (heat at open door) 4220 - 8440 


‘ Adapted from Table C-2-2.2.1(b) of Ref. 3 
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TABLE 3. RTI VALUES FOR ANY LISTED DETECTOR™ 
aeeee SS ee ee eee 


(FIXED TEMPERATURE TYPE 1/72 








Listed UL Listed Activation Temperature mit 
Spacing Listed 
(m) 53.3°C 57.2°C 62.8°C 71.1°C 76.7°C 91.1°C —_Temp. 
3.05 494 408 324 241 198 120 241 
4.57 309 235 193 136 110 56 136 
6.10 204 167 130 87 64 21 87 
7.63 153 124 96 59 40 59 
9.15 117 99 75 45 27 45 
12.20 88 71 51 7402 
S725 73 54 37 
2135 45 30 11 
NOTE: These RTIs are based on an analysis of the Underwriters 


Laboratories and Factory Mutual listing test procedures. Plunge test (see 
Ref. 7) results performed on the detector to be used will give a more 
accurate response time index. 


Abbreviations: 


eM 
UY 


Factory Mutual 
Underwriters Laboratories Inc. 


ue Adapted from Table C-3-2.1.1 of Ref. 5 
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TABLE 4. RTI VALUES FOR ANY LISTED DETECTOR” 
wee 
(RATE OF RISE TYPE) [(m-s)?/7] 


Listed UL Listed Activation Rate of Temp. Rise 
Spacing 
(m) 8.33°C/min 11.1°C/min 13.9°C/min 
3.05 1013 722 543 
3.81 802 593 445 
4.57 654 482 352 
6.10 482 321 235 
9.15 309 210 155 
12.20 247 161 114 
15.25 235 136 89 
Abbreviation: 
WG Underwriters Laboratories Inc. 


- Adapted from Ref. 3 


29 


S@AIN) YIMOIN PATZ 77 JO IOTd “T eain3t a 


(S) NOILINS! WOuS SWIL 
OOO! O06 O08 OOZL O09 OOS OOF OOE ONCE OOD! 





0 


(Mx) SLVY ASV3A1SY LVAH 


30 


uopAetez109 AZPOOTSA SseTUCFSUeT JO IOTd °Z eanstg 


H/4 





LO 
LO 


e1v//én 


Ol 


31 


This page intentionally left blank. 


32 


APPENDIX A. DETACT-T2 
FORTRAN VERSION 
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agaagaagagagggannnnagaagannngaaaannagaaagagagaaqaganaaaaannaa Q 


PROGRAM DETACT 


THIS IS PROGRAM DETACT-T2. 

IT CALCULATES THE RESPONSE OF THERMAL DETECTORS, BOTH FIXED TEMPERATURE 
AND RATE OF TEMPERATURE RISE, TO FIRES WHOSE HEAT RELEASE RATE GROWS 
PROPORTIONALLY WITH THE SQUARE OF TIME FROM IGNITION. 


WRITTEN BY DAVID W. STROUP, CENTER FOR FIRE RESEARCH, NBS 


THIS PROGRAM IS A CONTRIBUTION OF THE NATIONAL BUREAU OF STANDARDS AND 
IS NOT SUBJECT TO COPYRIGHT. 


VERSION 1.1 


RI AAR AR RRR AREAK NOTATION xi eseseter tek KKK AAA AAA 


INPUT: 
J 


UNITS CODE (1 OR 2) 

1 - INPUT DATA IN ENGLISH UNITS 
2 - INPUT DATA IN METRIC UNITS 
AMBIENT TEMPERATURE 

DETECTOR RESPONSE TIME INDEX 
DETECTOR ACTIVATION TEMPERATURE 
DETECTOR ACTIVATION RATE OF RISE 
CEILING HEIGHT 

DETECTOR SPACING 

GROWTH FACTOR CODE, CHARACTER VARIABLE (S, M, F, U, OR 0) 
FIRE GROWTH RATE CONSTANT 


TIME OF ACTIVATION FOR FIXED TEMPERATURE DETECTOR 
HEAT RELEASE RATE AT TIME OF ACTIVATION, T 

TIME OF ACTIVATION FOR RATE OF RISE DETECTOR 
HEAT RELEASE RATE AT TIME OF ACTIVATION, TR 
ERROR CODE (0 OR 1) 

0 - SUCCESSFUL 

1 - UNSUCCESSFUL 


KKK KAKAAKAKAAKAAAAKAKEKKAKRRARAREAR ARERR RRKRARERERRRERERERER 


IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
CHARACTER M 

DATA IRTTY, IWTTY/5,6/ 

DATA GE,CPE, RHOE/32.2D0,0.24D0,0.0735D0/ 
DATA GM,CPM,RHOM/9.8D0,1.0035D0,1.1768D0/ 
TERR = 0 

WRITE (IWTTY,1) 

FORMAT (' DETACT-T2 VERSION 1.0'/ 


& 
& 


' WRITTEN BY D.W. STROUP (1985)'// 
' CONTRIBUTION OF THE '/ 
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C 
C 
C 
5 


10 


20 


30 


40 


50 
60 
70 


80 


90 


100 


NATIONAL BUREAU OF STANDARDS (U.S.).'/ 
NOT SUBJECT TO COPYRIGHT. '// 
CALCULATES DETECTOR ACTUATION TIME'/ 
BELOW UNCONFINED CEILINGS WITH '/ 
TIME SQUARED FIRE GROWTH RATES.'//) 


R? BR | &S KS 


ENTER INPUT DATA 


CONTINUE 

WRITE (IWTTY,10) 

FORMAT (' ENTER: 1 FOR ENGLISH UNIT INPUT'/8X,'2 FOR METRIC UNIT I 
&NPUT'/) 

READ (IRTTY,*) J 

IF (J.EQ.1) THEN 

WRITE (IWTTY, 20) 

FORMAT (' ENTER THE AMBIENT TEMPERATURE IN DEGREES F.'/) 

READ (IRTTY,*) TAMB 

WRITE (IWTTY, 30) 

FORMAT (' ENTER THE DETECTOR RESPONSE TIME INDEX (RTI) IN (FT-SEC) 
&**1/2.'/) 

READ (IRTTY,*) RTI 

IF (RTI.LE.0.000001D0) RTI = 0.000001D0 

RTIR = RTI 

WRITE (IWTTY,40) 

FORMAT (' ENTER THE DETECTOR ACTIVATION TEMPERATURE IN DEGREES F.' 
&/) 

READ (IRTTY,*) TACT 

WRITE (IWTTY,50) 

FORMAT (' ENTER A DETECTOR RATE OF RISE IN DEGREES F/MINUTE. '/) 
READ (IRTTY,*) ROR 

WRITE (IWTTY,60) 

FORMAT (' ENTER THE CEILING HEIGHT IN FEET.'/) 

READ (IRTTY,*) HF 

WRITE (IWTTY,70) 

FORMAT (' ENTER THE DETECTOR SPACING IN FEET.'/) 

READ (IRTTY,*) ZF 

WRITE (IWTTY,80) 

FORMAT G ENTER: S FOR SLOW FIRE GROWTH RATE '/ 
& M FOR MEDIUM FIRE GROWTH RATE'/ 
& : F FOR FAST FIRE GROWTH RATE'/ 
& ‘ U FOR ULTRAFAST FIRE GROWTH RATE OR'/ 
& ' O FOR OTHER'/ ) 

READ (IRTTY,410) M 


IF ((M.EQ.'s').OR.(M.EQ.'S')) ALPHA = 0.00277778D0 

IF ((M.EQ.'m').OR.(M.EQ.'M')) ALPHA = 0.01111111D0 

IF ((M.EQ.'£').OR.(M.EQ.'F')) ALPHA = 0.0444445D0 

IF ((M.EQ.'u').OR.(M.EQ.'U')) ALPHA = 0.1777778D0 

IF ((M.EQ.'o').OR.(M.EQ.'O')) WRITE (IWTTY,90) 

IF ((M.EQ.'0').OR.(M.EQ.'O')) READ (IRTTY,*) ALPHA 

FORMAT (' ENTER THE FIRE GROWTH RATE CONSTANT (ALPHA) IN BTU/SEC/S 
&EC/SEC.'/) 

CONTINUE 
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C 


ELSE 


C METRIC UNIT INPUT 


C 


110 


130 


140 


150 


160 


170 


180 


C 


WRITE (IWTTY, 110) 

FORMAT (' ENTER THE AMBIENT TEMPERATURE IN DEGREES C.'/) 

READ (IRTTY,*) TAMB 

WRITE (IWTTY, 120) 

FORMAT (' ENTER THE DETECTOR RESPONSE TIME INDEX (RTI) IN (M-SEC)* 
&*1/2.'/) 

READ (IRTTY,*) RTI 

IF (RTI.LE.0.000001D0) RTI = 0.000001D0 

RTIR = RTI 

WRITE (IWTTY, 130) 

FORMAT (' ENTER THE DETECTOR ACTIVATION TEMPERATURE IN DEGREES C.' 
&/) 

READ (IRTTY,*) TACT 

WRITE (IWTTY,140) 

FORMAT (' ENTER A DETECTOR RATE OF RISE IN DEGREES C/MINUTE.'/) 
READ (IRTTY,*) ROR 

WRITE (IWTTY, 150) 

FORMAT (' ENTER THE CEILING HEIGHT IN METERS. '/) 

READ (IRTTY,*) HF 

WRITE (IWTTY, 160) 

FORMAT (' ENTER THE DETECTOR SPACING IN METERS.'/) 

READ (IRTTY,*) ZF 

WRITE (IWTTY,80) 

READ (IRTTY,410) M 


IF ((M.EQ.'s').OR.(M.EQ.'S')) ALPHA = 2.930555556D0 


IF ((M.EQ.'m').OR.(M.EQ.'M')) ALPHA = 11.72222222D0 

IF ((M.EQ.'£').OR.(M.EQ.'F')) ALPHA = 46.88888889D0 

IF ((M.EQ.'u').OR.(M.EQ.'U')) ALPHA = 187.5555556D0 

IF ((M.EQ.'o').OR.(M.EQ.'O')) WRITE (IWTTY,170) 

IF ((M.EQ.'o').OR.(M.EQ.'O')) READ (IRTTY,*) ALPHA 

FORMAT (' ENTER THE FIRE GROWTH RATE CONSTANT (ALPHA) IN WATTS/SEC/ 
&SEC.'/) 

CONTINUE 

ENDIF 


C weeded KA AKT CALCULATIONS xt reke dered de reset Rese KARIERRE IRE 


C 


R = 0.5D0*DSQRT(2.0D0)*ZF 
ROH = R/HF 

ROR = ROR/60.D0 

IF (J.EQ.1) THEN 

TAMB = TAMB + 460.D0 

TACT = TACT + 460.D0 

A = GE/(CPE*TAMB*RHOE) 
Gr=nGE 

ELSE 
TAMB 
TACT 


TAMB + 273.D0 
TACT + 273.D0 
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A = GM/(CPM*TAMB*RHOM*1000.D0) 
G = GM 
ENDIF 


TOTS2 = A**(-1.D0/5.D0)*ALPHA**(-1.D0/5.D0)*HF**(4.D0/5.D0) 
DLTODL = A**(2.D0/5.D0)*(TAMB/G)*ALPHA**(2.D0/5.D0) 


& *HF**(-3.D0/5.D0) 
UOUS2 = A**(1.D0/5.D0)*ALPHA**(1.D0/5.D0)*HF**(1.D0/5.D0) 
DELTD = TACT - TAMB 


IF (ROH.GT.0.3D0) THEN 

UODLTH = 0.59D0*ROH**(-0.63D0) 

ELSE 

UODLTH = 3.87D0/(9.115D0**0.5D0) 

ENDIF 

TS2F = 0.954D0*(1.D0+ROH) 

A2 = (4.D0/3.D0)*DLTODL*UOUS2**(-0.5D0)*UODLTH**(-0.5D0)*RTI/ 
& (TOTS2*(0.188D0+0.313D0*ROH) ) 

C = 1.DO+DELTD/A2 

CALL NWIN (C,Y,IERR, IWTTY) 

IF (IERR.EQ.1) GO TO 380 

DELTS2 = (4.D0/3.D0)*U0US2**(-0.5D0)*UODLTH**(-0.5D0)*RTI*Y/ 


& (TOTS2*(0.188D0+0.313D0*ROH) ) 
TS2 = 0.954D0*(1.D0+ROH)+(0.188D0+0. 313D0*ROH)*DELTS2**(3.D0/4.D0) 
TS2 = TS2 + TS2F 


T = TS2*A**(-1.D0/5.D0)*ALPHA**(-1.D0/5.D0)*HF**(4.D0/5.D0) 
QD = ALPHA*T**2 

IF (RTIR.LE.2.0D0) RTIR = 2.0D0 

D1 = (4.D0/3.D0)*DLTODL/ (TOTS2*(0.188D0+0.313D0*ROH) ) 

D2 = (3.D0/4.D0)*U0US2**0 . 5D0*UODLTH**0 . 5D0*(1.D0/RTIR)*TOTS2 
& *(0.188D0+0.313D0*ROH) 

CALL BISECT (D1,D2,ROR,DELTS2, IERR, IWTTY) 

IF (IERR.EQ.1) GO TO 380 

TS2R = 0.954D0*(1.D0+ROH)+(0.188D0+0.313D0*ROH)*DELTS2 
& **(3.D0/4.D0) 

TS2R = TS2R + TS2F 

TR = TS2R*A**(-1.D0/5.D0)*ALPHA**(-1.D0/5.D0)*HF**(4.D0/5.D0) 
QDR = ALPHA*TR**2 


C Keen PRINT OUT RESULTS Seve dete sede se cse seve sete sete fe se seve sete sete Fee 
C 

ROR = ROR*60.D0 

WRITE (IWTTY,190) 
190 FORMAT ('1l RESULTS:') 

IF (J.EQ.1) THEN 


C 
C ENGLISH (METRIC) OUTPUT 
Cc 
HF2 = HF * 0.3048D0 
ZF2 = ZF * 0.3048D0 


WRITE (IWTTY,200) HF,HF2,ZF,ZF2 

200 FORMAT ('0 CEILING HEIGHT ‘= ',F6.2,' FEET (',F6.2,' METERS)'/ 
& 5X, 'DETECTOR SPACING = ',F6.2,' FEET (',F6.2,' METERS)') 
RTI2 = RTI * (0.3048D0**0.5D0) 
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WRITE (IWTTY,210) RTI,RTI2 
210 FORMAT ('0 DETECTOR RTI = ',F6.1,' (FT-SEC)**1/2 (',F6.1,' (M-S 
&EC)**1/2)') 
ALPHA2 = ALPHA * 1055.D0 
WRITE (IWTTY,220) ALPHA, ALPHA2 
220 FORMAT ('0 FIRE GROWTH CONSTANT = ',E12.4,' BTU/(SEC**3)'/ 
& ' (',E12.4,' JOULES/SEC**3)') 
WRITE (IWTTY,400) 
READ (IRTTY,410) M 
WRITE (IWTTY, 230) 
230 FORMAT (/'O FOR TEMPERATURE ACTUATED DETECTOR: ' ) 
TACT = TACT - 460.D0 
TACT2 = (5.D0/9.D0)*(TACT-32.D0) 
WRITE (IWTTY,240) TACT, TACT2 
240 FORMAT ('0 ACTIVATION TEMPERATURE = ',F6.1,' DEGREES F (' ,F6.1, 
& ' DEGREES C)') 
T= T / 60.D0 
WRITE (IWTTY,250) T 
250 FORMAT ('0 TIME TO ACTIVATION = ',F8.2,' MINUTES') 
QD2 = QD * 1.055D0 
WRITE (IWTTY,260) QD,QD2 
260 FORMAT ('0 HEAT RELEASE RATE = ',E12.4,' BTU/SEC'/ 
& : (',E12.4,' KILOJOULES/SEC)') 
WRITE (IWTTY,400) 
READ (IRTTY,410) M 
WRITE (IWTTY, 270) 
270 FORMAT (//'O FOR RATE OF RISE ACTUATED DETECTOR: ' ) 
ROR2 = ROR * (5.D0/9.D0) 
WRITE (IWTTY,290) ROR,ROR2 
290 FORMAT ('0 ACTIVATION RATE OF RISE = ',F6.2,' DEGREES F/MIN (', 
& F6.2,' DEGREES C/MIN)') 
TR = TR / 60.D0 
WRITE (IWTTY,250) TR 
QDR2 = QDR * 1.055D0 
WRITE (IWTTY,260) QDR,QDR2 
WRITE (IWTTY,400) 
READ (IRTTY,410) M 


ELSE 
C 
C METRIC (ENGLISH) OUTPUT 
C 
HF2 = HF * (1.D0/0.3048D0) 
ZF2 = ZF * (1.D0/0.3048D0) 


WRITE (IWTTY,300) HF,HF2,ZF,ZF2 

300 FORMAT ('0 CEILING HEIGHT = ',F6.2,' METERS (',F6.2,' FEET)'/ 
& 5X, 'DETECTOR SPACING = ',F6.2,' METERS (',F6.2,' FEET)') 
RTI2 = RTI * ((1.D0/0.3048D0)**0.5D0) 
WRITE (IWTTY,320) RTI,RTI2 

320 FORMAT ('0 DETECTOR RTI = ',F6.1,' (M-SEC)**1/2 (',F6.1,' (FT-S 
&EC)**1/2)') 
ALPHA2 = ALPHA * (1.D0/1055.D0) 
WRITE (IWTTY,330) ALPHA,ALPHA2 
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330 


340 


350 


360 


370 


C 


Cc 


FORMAT (' FIRE GROWTH CONSTANT = ',E12.4,' JOULES/(SEC**3)'/ 
& ' (',E12.4,' BTU/SEC**3)') 
WRITE (IWTTY,400) 

READ (IRTTY,410) M 

WRITE (IWTTY, 230) 

TACT = TACT - 273.D0 

TACT2 = (9.00/5.D0) * TACT + 32.D0 

WRITE (IWTTY,340) TACT, TACT2 

FORMAT ('0 ACTIVATION TEMPERATURE = ',F6.1,' DEGREES C (',F6.1, 
& ' DEGREES F)') 

T=T / 60.D0 

WRITE (IWTTY,350) T 

FORMAT ('0 TIME TO ACTIVATION = ',F8.2,' MINUTES') 

QD = QD / 1000.D0 

QD2 = QD * (1.D0/1.055D0) 

WRITE (IWTTY, 360) QD,QD2 
FORMAT ('0 HEAT RELEASE RATE = ',E12.4,' KILOJOULES/SEC'/ 
& ‘ (',E12.4,' BTU/SEC)') 
WRITE (IWTTY, 400) 

READ (IRTTY,410) M 

WRITE (IWTTY, 270) 

ROR2 = ROR * (9.D0/5.D0) 

WRITE (IWTTY,370) ROR,ROR2 

FORMAT ('0 ACTIVATION RATE OF RISE = ',F6.2,' DEGREES C/MIN (', 
& F6.2,' DEGREES F/MIN)') 

TR = TR / 60.D0 

WRITE (IWTTY,350) TR 

QDR = QDR / 1000.D0 

QDR2 = QDR * (1.D0/1.055D0) 

WRITE (IWTTY,360) QDR,QDR2 

WRITE (IWTTY, 400) 

READ (IRTTY,410) M 

ENDIF 

WRITE (IWTTY, 420) 

READ (IRTTY,410) M 

IF ((M.EQ.'Y').OR.(M.EQ.'y')) GO TO 5 

STOP 'PROGRAM COMPLETED' 

CONTINUE 

WRITE (IWTTY, 390) 

FORMAT (' *%%*#%%%> ERROR IN DETACT ROUTINE <**#%#%#%!) 

FORMAT (//' <RETURN> TO CONTINUE' ) 

FORMAT (Al) 

WRITE (IWTTY,420) 

FORMAT (' TRY ANOTHER CASE (Y/N) ?'/) 

READ (IRTTY,410) M 

IF ((M.EQ.'Y').OR.(M.EQ.'y')) GO TO 5 

STOP 'PROGRAM ABORTED' 

END 


SUBROUTINE NWIN (C,P,IERR, IWTTY) 


C THIS IS ROUTINE NWIN. IT SOLVES THE FIXED TEMPERATURE DETECTOR 
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C ACTIVATION TIME EQUATION USING A NEWTON-RAPHSON METHOD. 
C 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
PO = 0.1D0 
TOL = 0.00001D0 
NO = 1000 
r= 1 
IERR = 0 
10 CONTINUE 
IF (I1.LE.NO) THEN 
IF (PO.GT.50.D0) THEN 
X = 0.0D0 
ELSE 
X = DEXP(-PO) 
ENDIF 
FX = PO+X-C 
FPMX = 1.0D0 - X 
IF (FPMX.LT.0.00000001D0) GO TO 30 
P = PO - (FX/FPMX) 
IF (ABS(P-P0).LT.TOL) THEN 
IERR = 0 
RETURN 
ELSE 
Sa <p i 
PO = P 
ENDIF 
GO TO 10 
ENDIF 
TERR = 1 
WRITE (IWTTY,20) I 
20 FORMAT (' NEWTON-RAPHSON FAILED AFTER NO ITERATIONS, NO = ',I4) 
RETURN 
30 WRITE (IWTTY, 40) 
40 FORMAT (' SLOPE OF EQUATION TOO CLOSE TO ZERO FOR '/ 


& '  NEWTON-RAPHSON METHOD. '/ 
& ' ERROR RETURN' ) 
TERR = 1 
RETURN 
END 
C 
SUBROUTINE BISECT (D1,D2,ROR,P, IERR, IWITY) 
C 
C THIS IS ROUTINE BISECT. IT SOLVES THE RATE OF RISE DETECTOR ACTIVATION 
C TIME EQUATION USING A BISECTION METHOD. 
C 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
IERR = 0 
TOL = 0.00001D0 
NO = 1000 
A = 0.0D0 
B = 1000.0D0 


RLMT = TOL / 20.0 
40 


10 


20 


30 


CONTINUE 
IF ((D2*B).GT.50.D0) THEN 
X = 0.0D0 
ELSE 
X = DEXP (-D2*B) 
ENDIF 
FXB = D1*B**0.25D0-D1*B**0.25*X-ROR 
IF (FXB.LT.0.0D0) THEN 
B 
B + 500.D0 
0 TO 10 

ENDIF 

i 

CONTINUE 

IF (I1.LE.NO) THEN 

P = A+(B-A)/2.0D0 

IF ((D2*P).GT.50.D0) THEN 
Xe=OOD0 

ELSE 

X = DEXP (-D2*P) 

ENDIF 

FX = D1*P**0.25D0-D1*P**0.25D0*X-ROR 
IF (((FX.GT.-RLMT).AND.(FX.LT.RLMT)).OR. 
&  (((B-A)/2.D0).LT.TOL)) THEN 
TERR = 0 

RETURN 

ELSE 

Tes ye 4h i 

IF ((D2*A).GT.50.D0) THEN 

a= ORONO 

ELSE 

X = DEXP (-D2*A) 

ENDIF 

FXA = D1*A**0.25D0-D1*A**0.25D0*X-ROR 
IF ((FXA*FX).GT.0.0D0) THEN 
A=P 

ELSE 

B=P 

ENDIF 

ENDIF 

GO TO 20 

ENDIF 

IERR = 1 

WRITE (IWTTY,30) I 


A 
B 
G 


FORMAT (' BISECTION ROUTINE FAILED AFTER NO ITERATIONS, NO = ',14) 


RETURN 
END 
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APPENDIX B. DETACT-T2 
BASIC VERSION 


43 


PRINT "DETACT-T2 VERSION 2.0" 
PRINT "WRITTEN BY D.W. STROUP (1985)" 
PRINT "' "' 
PRINT "CONTRIBUTION OF THE" 
PRINT "NATIONAL BUREAU OF STANDARDS (U.S.)." 
PRINT "NOT SUBJECT TO COPYRIGHT. " 
PRINT "' "' 
PRINT "CALCULATES DETECTOR ACTUATION TIME" 
PRINT "BELOW UNCONFINED CEILINGS WITH " 

PRINT "TIME SQUARED FIRE GROWTH RATES " 

PRINT '" "' 

REM 

REM INPUT 

REM J - UNITS CODE (1 OR 2) 

REM 1 - INPUT DATA IN ENGLISH UNITS 

REM 2 - INPUT DATA IN METRIC UNITS 

REM TAMB - AMBIENT TEMPERATURE 

REM RTI - DETECTOR RESPONSE TIME INDEX 

REM TACT - DETECTOR ACTIVATION TEMPERATURE 

REM ROR’ - DETECTOR ACTIVATION RATE OF RISE 

REM HF - CEILING HEIGHT 

REM ZF - DETECTOR SPACING 

REM M - GROWTH RATE FACTOR CODE, CHARACTER (S, M, F, U, OR 0) 
REM ALPHA - FIRE GROWTH RATE CONSTANT (FOR M = O) 

REM 

REM OUTPUT: 

REM T - TIME OF ACTIVATION FOR FIXED TEMPERATURE DETECTOR 
REM QD - HEAT RELEASE RATE AT TIME OF ACTIVATION, T 
REM TR - TIME OF ACTIVATION FOR RATE OF RISE DETECTOR 
REM QDR  ~- HEAT RELEASE RATE AT TIME OF ACTIVATION, TR 
REM IERR - ERROR CODE (0 OR 1) 

REM 0 - SUCCESSFUL 

REM 1 - UNSUCCESSFUL 

REM 

CLEAR 

REM 

REM ENTER INPUT DATA 

REM 

READ GE#,CPE# , RHOE}# 

DATA 32.2#,0.24#,0.0735# 

READ GM#,CPMi# , RHOM# 

DATA 9.8#,1.0035#,1.1768# 

TERRZ = 0 

PRINT '"' ENTER: 1 FOR ENGLISH UNIT INPUT" 

PRINT "' 2 FOR METRIC UNIT INPUT" 

INPUT "SELECTION = ",JZ 

IF J%<>1 GOTO 740 

REM 

REM ENTER DATA IN ENGLISH UNITS 

REM 
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= 


510 INPUT "AMBIENT TEMPERATURE <DEGREES F> = ", TAMB# 

520 INPUT "DETECTOR RESPONSE TIME INDEX, RTI <(FT-SEC) (1/2)> = ",RTI# 
530 IF RTI#<.0000001/ THEN LET RTI# = .0000001} 

540 RTIR# = RTI# 

550 INPUT "DETECTOR ACTIVATION TEMPERATURE <DEGREES F> = ", TACT 

560 INPUT "DETECTOR RATE OF TEMPERATURE RISE <DEGREES F/MIN> = ",ROR# 
570 INPUT "ROOM CEILING HEIGHT <FEET> = ",HF# 

580 INPUT "DETECTOR SPACING <FEET> = ",,ZF# 


590 PRINT "ENTER: S FOR SLOW FIRE GROWTH RATE," 
600 PRINT " M FOR MEDIUM FIRE GROWTH RATE," 
610 PRINT "' F FOR FAST FIRE GROWTH RATE," 
620 PRINT "' U FOR ULTRAFAST FIRE GROWTH RATE, OR" 
630 PRINT "' O FOR OTHER FIRE GROWTH RATE" 


640 INPUT "SELECTION = ",M$ 


650 IF (M$="s" OR M$="S") THEN LET ALPHA# = .00277778# 
660 IF (M$="m" OR M$="M'"') THEN LET ALPHA# = .011111111 
670 IF (M$="£" OR M$="F'') THEN LET ALPHA# = .04444445} 
680 IF (M$="u"" OR M$="U'"') THEN LET ALPHA# = .17777778}# 


690 IF (M$="0" OR M$="0"') THEN INPUT "FIRE GROWTH RATE <BTU/SEC 3> = ",ALPHA# 
700 GOTO 940 

710 REM 

720 REM ENTER DATA IN METRIC UNITS 

730 REM 

740 INPUT "AMBIENT TEMPERATURE <DEGREES C> = ", TAMB# 

750 INPUT "DETECTOR RESPONSE TIME INDEX, RTI <(M-SEC) (1/2)> = ",RTI# 
760 IF RTI##<.00000014 THEN LET RTI# = .0000001# 

770 RTIR# = RTI# 

780 INPUT "DETECTOR ACTIVATION TEMPERATURE <DEGREES C> = ", TACT? 

790 INPUT "DETECTOR RATE OF TEMPERATURE RISE <DEGREES C/MIN> = ",ROR# 
800 INPUT ''ROOM CEILING HEIGHT <METERS> = ",HF# 

810 INPUT "DETECTOR SPACING <METERS> = ", ZFif 


820 PRINT "ENTER: S FOR SLOW FIRE GROWTH RATE," 
830 PRINT "' M FOR MEDIUM FIRE GROWTH RATE," 
840 PRINT "' F FOR FAST FIRE GROWTH RATE," 
850 PRINT "' U FOR ULTRAFAST FIRE GROWTH RATE, OR" 
860 PRINT "' O FOR OTHER FIRE GROWTH RATE" 


870 INPUT ''SELECTION = ",M$ 

880 IF (M$="s" OR M$="S") THEN LET.ALPHA# 

890 IF (M$="m" OR M$="M") THEN LET ALPHA# 

900 IF (M$="£" OR M$="F") THEN LET ALPHA# = 46.888889} 

910 IF (M$="u" OR M$="U") THEN LET ALPHA# = 187.55556# 

920 IF (M$="0" OR M$="0") THEN INPUT "FIRE GROWTH RATE <J/SEC 3> = ",ALPHA# 
R 


2.9305556# 
11.722222# 


QAO REM xed ~~ CALCULATIONS xx eve redesedetede ede tee tee HAHAHAHAHA RRR 


950 REM 
960 R# = .5i# * SQR(2#) * ZF 

Ri /HF# 

ROR# / 60# 

990 IF J%<>1 GOTO 1050 

1000 TAMBi# = TAMB# + 460% 

1010 TACT# = TACT# + 460% 

1020 Att = GE#/(CPE#*TAMB#*RHOE} ) 


V2) 
co 
j=) 
2] 
[e} 
be) 
= 
oul 
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Git = GE# 

GOTO 1090 

TAMB/## = TAMB# + 273# 
TACT# = TACT# + 273+ 
Aff = GM#/ (CPM##*TAMB#*RHOM#* 10004) 
Gi# = GM# 

TOTS2#=A# (-14/54)*ALPHA# (-14/54)*HFH (44/54) 

DLTODL#=A# (2d / 54) * (TAMB#/ Git) *ALPHA# (2d / 5d) *HF# (-34/ 544) 
UOUS2#=Ad (14/5) *ALPHAH (14/54)*HFH# (14/54) 
DELTD# = TACT# - TAMBi# 
IF (ROH#>.3) THEN UODLTH#=.59*ROH# (-.63/) ELSE UODLTH#=3.87#/(SQR(9.115#) ) 

TS2F# = .954i * (14 + ROH#) 

PH = TOTS2#*(. 188+. 3134*ROH#) 

A2s#H=(4#/ 34) *DLTODL#*UOUS2# (- . 54)*UODLTH# (- . 5##) *RTI#/ Pit 

CH = 14 + DELTD# / A2# 

GOSUB 2650 

IF IERR%=1 GOTO 2550 

DELTS2#=( 4#/3#)*UOUS2#° (- . S#)*UODLTH# (- . 54) *RTI#*Y##/ Pit 

TS2#=. 544% (1#+ROH#H)+(. 188+. 3134*ROH#)*DELTS2#° (34 / 4#) 

TS2# = TS2# + TS2F# 

TH=TS2#* A (- 1/54) *ALPHAH (-14t/ Sit) *HEH (44/ 5H) 

QD/t = ALPHA# * Ti# 2! 

IF RTIR#<=2i# THEN LET RTIR# = 2+ 

D1f=(4#/ 34) *DLTODL# / Pit 

D2#=( 34/44) *UOUS2# . S#UODLTH# . 54##*( 14 /RTIR#) * PH 

GOSUB 2950 

IF IERRZ=1 GOTO 2550 

TS2R#t=. 9544 ( 144+ROH#)+( . 1884+. 313i#*ROH#)*DELTS2# (34 / 4H) 

TS2R# = TS2R# + TS2F# 

TRH=TS2R#H*AH (- 14/54) *ALPHAF (- 14/54) *HFH (4d / 5i#) 








QDR# = ALPHA# * TR# 2! 

ROR# = ROR# * 60# 

REM 

REM PRINT OUT RESULTS 

REM | 

PRINT "RESULTS:" 

IF J%<>1 GOTO 1960 

REM 

REM ENGLISH (METRIC) OUTPUT 

REM 

HF2# = HF# * .3048 

ZF2# = ZF# * .3048 

PRINT "' "' 

PRINT '' CEILING HEIGHT = ';CSNG(HF#);"' FEET" 
PRINT " ('"; CSNG(HF2#) ;'' METERS )" 
PRINT "' " 

PRINT " DETECTOR SPACING = "';CSNG(ZF#);"' FEET" 
PRINT "' ('"; CSNG(ZF2#) ;"' METERS )" 
PRINT" "' 

RTI2# = RTI# * (.3048#° . 5H) 


PRINT " RTI = ";CSNG(RTI#);" (FT-SEC) 0.55" ('"3CSNG(RTI2#);" (M-SEC) 0.5) 
PRUN Tee 
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ALPHA2/ = ALPHA# * 1055# 


PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 
PRINT 
TACT#H 


it 


i] 


FIRE GROWTH RATE = ";CSNG( ALPHA?) ;"" BTU/(SEC 3)" 
(''; CSNG(ALPHA2#) ;"" J/(SEC 3))" 
tl 
<RETURN> TO CONTINUE ";0$ 
in} 


FOR TEMPERATURE ACTUATED DETECTOR:" 
TACT# - 460% 


TACT2# = (5#/9#) * (TACTH# - 32+) 


PRINT 
PRINT 
PRINT 
PRINT 


ACTIVATION TEMPERATURE = ";CSNG(TACT#);" DEGREES F" 
("";CSNG(TACT2#);" DEGREES C)" 


TH = Ti / 604 
PRINT " 
PRINT AL! 

QD2# = QD# * 1.0554 


PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 
PRINT 
ROR2# 
PRINT 
PRINT 
PRINT 
PRINT 
TR# = 
PRINT 
QDR2# 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


T 


at 


TIME TO ACTIVATION = '';CSNG(Ti#);" MINUTES" 


HEAT RELEASE RATE = ";CSNG(QD#);' BTU/SEC" 
(''s CSNG(QD2#);'' kJ/SEC )" 
in} 
<RETURN> TO CONTINUE "';0$ 
it 


FOR RATE OF RISE ACTUATED DETECTOR:" 
ROR# * (5#/9#) 
i] 


ACTIVATION RATE OF RISE = ";CSNG(ROR#);'' DEGREES F/MIN" 


('"; CSNG(ROR2#);'' DEGREES C/MIN)" 


Rit / 60+ 
TIME TO ACTIVATION = ";CSNG(TR#);"' MINUTES" 
QDR# * 1.055# 


HEAT RELEASE RATE = ";CSNG(QDR#);"' BTU/SEC" 
(''; CSNG(QDR2#);'' kJ/SEC )" 


it 


INPUT '' TRY ANOTHER RUN (YES/NO) '';0$ 
IF (O$="Y" OR O$="'y"') GOTO 120 
2530 


GOTO 


HF# * (1!/.3048) 
ZF * (1!/.3048) 


REM METRIC (ENGLISH) OUTPUT 


" 


CEILING HEIGHT = '';CSNG(HF#);'' METERS" 
(''; CSNG(HF2#) ;"" FEET)" 
10 
DETECTOR SPACING = ";CSNG(ZF#);'' METERS" 
(''s CSNG(ZF2#);'' FEET)" 
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PRINT "' "' 

RTI2# = RTI# * ((1!/.3048)° .5#) 

PRINT '"' RTI = '';CSNG(RTI#);"' (M-SEC)°0.5"';" ('"';CSNG(RTI2#) 3" (FT-SEC) 0.5) 
PRINT '"' "' 

ALPHA2/## = ALPHA# * (1!/1055!) 

PRINT "' FIRE GROWTH RATE = ";CSNG(ALPHA#);' J/(SEC 3)" 

PRINT " (''; CSNG(ALPHA2#);'' BTU/(SEC 3))" 
PRINT "' " 

INPUT '<RETURN> TO CONTINUE ";0$ 

PRINT "' " 

PRINT "' "' 

PRINT '' FOR TEMPERATURE ACTUATED DETECTOR:" 

PRINT "' " 

TACT# = TACT# - 273+ 

TACT2# = (9f/5#) * TACT# + 32¢# 

PRINT '' ACTIVATION TEMPERATURE = '";CSNG(TACT#);'' DEGREES C" 
PRINT "' (";CSNG(TACT2#);'' DEGREES F)" 
PRINT " " 

TH = Tit / 604 

PRINT '' TIME TO ACTIVATION = ''s;CSNG(T#);'' MINUTES" 

PRINT " "' 

QDi# = QD# / 1000# 

QD2# = QD# * (1!/1.055) 

PRINT '' HEAT RELEASE RATE = '"';CSNG(QD#);"' kJ/SEC" 


PRINT "' (''sCSNG(QD2#);'" BTU/SEC )" 
PRINT "' "' 

INPUT "<RETURN> TO CONTINUE"; 0$ 

PRINT "' "' 

PRINT '' '"! 

PRINT '' FOR RATE OF RISE ACTUATED DETECTOR: "' 


ROR2# = ROR# * (9#/5i#) 
PRINT "' "' 

PRINT ' ACTIVATION RATE OF RISE = ";CSNG(ROR#);'' DEGREES C/MIN" 
PRINT "' (''; CSNG(ROR2#);'' DEGREES F/MIN)" 
PRINT " "' 

TRH = TR# /60# 

PRINT ' TIME TO ACTIVATION = ''s;CSNG(TR#);'' MINUTES" 
PRINT " " 

QDR# = QDR# / 1000+ 

QDR2# = QDR# * (1!/1.055) 

PRINT '' HEAT RELEASE RATE = ''sCSNG(QDR#);"' kJ/SEC" 
PRINT ('';CSNG(QDR2#);'' BTU/SEC )" 
PRINT "' "' 

PRINT " " 

INPUT '' TRY ANOTHER RUN (YES/NO) '';0$ 

IF (O$="Y"" OR O$="y") GOTO 120 

END 

REM 

REM Sekekededetecesesesese ERROR HANDLING ROUTINE KRHA 

REM 

PRINT "' '"' 

PRINT "' !!!!1 ERROR IN DETACT-T2 ROUTINE !!!!!" 
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PRINT " !!!!! PROGRAM RUN ABORTED ue 
PRON 
INPUT '' WOULD YOU LIKE TO TRY AGAIN (YES/NO) ?'';0$ 


IF (O$="Y" OR O$="y") GOTO. 120 

END 

REM 

REM ** NEWTON RAPHSON ROUTINE TO FIND ACTIVATION ** 
REM ** TIME OF FIXED TEMPERATURE DETECTOR. we 
REM 

PO = .1# 

TOL# = .00001# 

NOZ = 1000 

IZ = 

IERRZ = 0 

REM 

IF IZ>NOZ GOTO 2870 


IF (POf#>50#) THEN LET Xf=0f ELSE X##=EXP(-PO#) 
FX#=PO# + Xi# - CH 

FPMX# = lit - Xi# 

IF (FPMX#<.00000001#) GOTO 2900 
P2# = PO# - (FXi#/FPMX#) 

IF (ABS(PO#-P2#)>TOL#) GOTO 2840 
IERRZ%Z = 0 

Y# = P2# 

RETURN 

1 2 164 & il 

PO# = P2# 

GOTO 2730 

ITERRZ = 


PRINT "NEWTON RAPHSON FAILED AFTER NO ITERATIONS, NO = 


RETURN 
PRINT "SLOPE OF EQUATION TOO CLOSE TO ZERO!!!" 
PRINT "ERROR IN NEWTON RAPHSON ROUTINE" 


IERRZ = 1 

RETURN 

REM 

REM ** BISECTION ROUTINE TO FIND ACTIVATION ** 
REM ** TIME OF RATE OF RISE DETECTOR. wy 
REM 

IERRZ = 

TOL# = .00001# 

NOZ = 1000 

AA# = Of 

Bi = 1000} 


RLMT# = TOL# / 20} 


REM 
IF (D2#*Bi#)>50/ THEN Xf = Off ELSE X# = EXP(- D2#*B#) 
FXB# = Diff * Bit -25# - D1dt * Bit 254 * Xft - ROR# 
IF FXBf>=0# GOTO 3110 

AA# = Bit 

Bit = Bit + SOO0# 

GOTO 3040 
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Mina 4 


3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 


TZ = 1 

REM 

IF I>NO GOTO 3260 

P3# = AA + (BH - AA) / 24 

IF (D2#*P3)>504# THEN X# = Off ELSE X# = EXP (-D2#*P3#) 
FX# = Diff * P3f ~ .25f - Dlé#t * P3H# ~ .25i * X# - ROR# 
DELTS2# = P3# 

IF ((FX#>(-RLMT#)) AND (FXi#<RLMT#)) THEN RETURN 

IF ((Bé-AA#) /2é#)<TOL# THEN RETURN 

REM 

IZ=I1Z%+1 

IF (D2/HAAd#)>50# THEN Xf = Of ELSE X## = EXP(-D2/*AA#) 
FXA# = Dli#f * AA# ~ .25¢ - Didf * AA# ~ .25¢ * X/# - ROR# 
IF (FXA#*FX#)>O# THEN AA# = P3/ ELSE Bi = P3# 

GOTO 3120 

IERRZ = 1 

PRINT "BISECTION ALGORITHM FAILED TO CONVERGE AFTER ",1Z%,' ITERATIONS." 
RETURN 
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DETECTOR ACTIVATION TIME TABLES: 
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